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Detection of continuum radio emission associated with Geminga 
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ABSTRACT 

A deep Very Large Array observation of the Geminga pulsar field led to the discovery, at a 
higher than lOcr significance level, of radio emission trailing the neutron star proper motion. 
This ~10''-long radio feature, detected with a flux of ~0.4 mJy at 4.8 GHz, is marginally 
displaced (2.7 ± 1.8 arcsec) from the pulsar (which, at any rate, is unlikely to contribute 
with magnetospheric pulsed emission more than 15% to the total observed radio luminosity, 
~10^^ erg s~^) and positionally coincident with the X-ray axial tail recently discovered by 
Chandra and ascribed to the pulsar wind nebula (PWN). Overall, the Geminga radio tail is 
compatible with the scenario of a synchrotron-emitting PWN, but the present data do not allow 
us to discriminate between different (and not always necessarily mutually exclusive) possible 
processes for producing that. If this radio feature does not result from intrinsic peculiarities 
of Geminga, but its proximity and radio-quiet nature (both helping in not hindering the faint 
diffuse radio emission), other relatively near and energetic radio-quiet pulsars could show 
similar structures in dedicated interferometric observations. 

Key words: pulsars: general - pulsars: individual (J0633+1746, Geminga) - radio contin- 
uum: stars - stars: neutron - stars: winds, outflows. 



1 INTRODUCTION 

Geminga (PSR J0633+1746) is a nearby radio-quiet pulsarQ dis- 
covered as a gamma-ray source and later identified as a neu- 
tron star (NS) th rough optical and X-ray observations (see 
iBignami & Caraveo 1996 for review). The pulsar period (P ^237 
ms) and its derivative (P ~1.1 x 10~^^ s s~^) corresponcQ to a 
spin-down age r ^^340 kyr, a spin-down power E^ot = 3.3 x 10^^ 
erg s~^ and a surface magnetic field Psurf ^1.6x10^^0. Hubble 
Space Telescope parallax measurements confirmed the proximity of 
Geminga {d — 250^52° P^) ^ pr oper motion correspo nding to 
a transverse velocity of ^210 km s~^ Jpahertv et al .l2007h , exceed- 
ing the typical sound speed in the interstellar medium (10-30 km 
s~^). Together with Crab and Vela, Geminga is one of known pul- 
sars with the highest spin-down flux {Erotd~^), allowing detailed 
studies of weak structures in the immediate vicinity of the NS and 
its interactions with the local interstellar medium. 

Pulsar wind nebulae (PWNe) are common diffuse features 
surrounding NSs with very different spin-down ages, that result 
from the interaction of the relativistic pulsar wind and the ambi- 
ent medium, producing shocks and outflows that can be observed, 
mostly because of synchrotron and inverse Compton (IC) radia- 
tion, in a very broad range from radio to gamma-ray energies ( see 



E-mail: apellizz@oa-cagliari.inaf.it 
^ We refer to a pulsar as radio-quiet when, even though it may emit radio 
pulses, these cannot be detected at Earth. 

^ See e.g. Lorimer & Kra mer 2005 for the formulae used to derive r, E^rot 
and ^surf from P and P. 



iKaspi et al.' "20 061: iGaensler & Stand l2006l : iKargaltsev & Pavlovl 
2008; Peflizzon ietal.ll201Qh . 

The complex PWN structures mostly depend on the partic- 
ular wind outflow geometry and on the r atio of the pulsar speed 
to the sound speed in the ambient medium feucciantini et al.ll2005l : 
iRomanova et al .l2005l : lBogovalov et al.l2005l\ Well-resolved equa- 
torial "tori" and axial (with respect to the spin axis) "jets" are 
seen around pulsars moving in the interste llar medium with sub- 
sonic velocities (e.g . Crab and Vela pulsars; IWeisskopf et"aDl200Ql : 
iHelfand et al.l200lh . while in pulsar moving with supersonic speed, 
as is the case of Geminga, bow-like shapes ahead of the pulsar and 
tails behind are typically observed (Pellizzoni et al. 2005 ). 

In fact, XMM-Newton observations of Geminga in April 2002 
revealed two ^2^ -long tails behind the pulsar, approximately sym- 
metri c with respect to th e sky projection of the pulsar's trajec- 
tory (ICaraveo et al.ll2003l) . Geminga was then observed by Chan- 
dra in February 2004 unveiling a new structure ^25^^ -long and 
^5^^-thick, with a surface brightness ^40 times higher than that of 
the XMM-Newton tails, starting at t he pul sar position and aligned 
with the proper motion direction (jPe Lu ca et al. 2006). A faint 
arc- Hke structure was also reported 5-7 arcsec ahead of the pul- 
sar jPavlov et al.ll2006h . A deeper Chandra observation carried out 
in August 2007 confirmed the existence of the axial tail (50^^- 
long) and of the two outer tails discovered by XMM-Newton 
JPavlov et al.ll201Ql) providing an overall X-ray PWN luminosity 
of Lo.3-8 kev = 3 X lO^^diso crg s-^ (^10_fog^ii^rot). Com- 
paring 2004 and 2007 Chandra observations, IPavlov et all (1201 Oh 
found indication that the X-ray tail is composed of variable ~5^^- 
long sub- structures (streaming "blobs"). This complex morphology 
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was interpreted as synchrotron radiation possibly arising from the 
shocked pulsar wind collimated by the ram pressure and/or jets em- 
anating along the pulsar's spin axis. 

Apart from diffuse TeV emission around Geminga on larger 
scales (>2°) with respect to the X-ray nebula dAbdo et al.ll2009h 
and a possible detection in the mid-infrared (jPanilenko et all 
l201lh . the Geminga PWN has not been detected in other energy 
bands so far. In particular, since the discovery of Geminga as a 
gamma-ray source, many observers have attempted to detect it at 
radio frequencies, both as a continuum and as a pulsating source, 
providing only tentative detection s of weak pulsed emission at 
low frequencies (^100 MHz: see [ Kuzmin & Losovskvl 1997allbl: 
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Malofeev & M alov 1997; Shitov et al. 1997; 'Shitov & PugachevI 
1998l : IVats et al.J997 . 1999) . Kassim & Lazio ( 1999) reported up- 



per limits on the continuum emission (Very Large Array [VLA] 
data) at 74 MHz «56 mJy) and 326 MHz «5 mJy) reconcil- 
ing their nondetection with previous pulsed detections by invoking 
intrinsic or extrinsic (refra ctive interstellar scinti l lation) variabil- 
ity. At higher frequencies, Ispoelstra & HermsenI (Il984h reported 
an upper limit of 0.5 mJy at 21 cm (Westerbork Synthesis Radio 
Telescope data) and 1 mJy at 6 cm (VLA data). More recently, 
iGiacani et al ](l2QQl) presented a high-resolution (24^^) study of the 
Hi interstellar gas distribution around Geminga including a 1.4 
GHz upper limit of ~0.4 mJy in the pulsar direction. Here we 
present the results of the deepest VLA interferometric observation 
of Geminga performed so far. 



2 VLA OBSERVATIONS AND DATA ANALYSIS 

Geminga was observed at 4.8 GHz (6 cm), with a bandwidth of 
50 MHz with the VLA in D configuration (program AP468). The 
observation was performed on 2004 July 24 for a total of 6.0 
on-source hours, with the antennas pointed at R.A.(J2000) = 
06^33™54!02, Dec.(J2000) = 17°46'll'.'5. Calibration and 
imaging were performed with the NRAO Astronomical Image 
Processing System (AlPS) software packaged The flux-density 
scale was calibrated by observing 0137-1-331 (3C48). The source 
0625+146 was observed at intervals of 30 minutes and used as 
phase calibrator. The surface brightness image was produced fol- 
lowing the standard procedures: Fourier-Transform, Clean, and Re- 
store implemented in the AIPS task IMAGR. We averaged the 2 IPs 
together in the gridding process under IMAGR. Self-calibration was 
applied to remove residual phase variations. 

In Fig. [T] we show a 8^ -side field of the radio contin- 
uum emission around Geminga. The image has a FWHM beam 
of 15'.'2 X 11^8 (position angle, PA = 46.5°) and a noise 
level of 25 /iJy beam~^ (la). A radio feature, with a peak of 
brightness at >10cr significance level, located at R.A.(J2000) = 
06^33^54^22 (±0^04), Dec.(J2000) = 17°46'll'.'22 (±0'.'42) 
(through AIPS task JMFIT), has been detected. The absolute astro- 
metric accuracy of VLA observations is linked to the uncertainty 
in the phase calibrator position, which is better than 0^.^15 in our 
case. However, the VLA accuracy can be limited by a number of 
effects, such as the atmospheric phase stabilityQFor this reason, we 
adopted for the positioning uncertainty a more conservative figure 
of 1^.^8, which was estimated as the synthesised beam size divided 
by the signal to noise ratio (see, e.g. jHagiwara et al1l200lh . 
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Figure 1. Geminga field (~8^ x 8^ at 4.8 GHz resulting from the 6- 
hours VLA observation performed in D configuration. The beam FWHM 
is 15(^2 X 11(^8 and the image sensitivity is 25 fiiy beam~^ (Icj). A zoom 
of the image, restored with a circular beam of 15.3^^ x 15.3^^ is inset in 
order to evidenciate the slightly elongated nature of the radio feature. The 
faintest drawn radio contour level corresponds to 75 //Jy beam~^ and the 
other levels are spaced by a factor of 



The observed brightness distribution is more extended than 
that expected for a point source: it is slightly elongated and fits with 
a single 2-dimensional Gaussian model with a FWHM major axis 
of 18V6±ir6,aminoraxisof liridzirO,andPA = (64.3db6.6)°, 
incompatible with beam parameters (the fit with the point- spread 
function yields a poor reduced of ^^2.5 for 5 degrees of free- 
donB. The deconvolved size along the major axis (approximately 
the direction of Geminga proper motion) is 11^^ =1= 3^^ while it 
is unresolved in the transverse direction. The flux density of the 
source at 4.8 GHz is (0.37 ± 0.05) mJy being calculated by in- 
tegrating the total intensity surface brightness down to the noise 
level. No evidence of source variability was found down to 1.5-2 
hours, which is the minimum time-scale that can be probed in our 
data. At the distance of Geminga, the corresponding radio power is 
^-SGHz - 2.7 X 10^ WHz-\ 

To investigate the presence of such source in other VLA ob- 
servations, we analyzed an archival dataset at 1.5 GHz in B config- 
uration (program AK147). In this case we obtained an image with 
a FWHM beam of 4^ 10 x 3^98 and noise level of 22 /iJy beam" ^ 
(Icr). We did n ot found any source in correspondence to the 4.8 
GHz detection. iGiacani et all (l2005h reported on a study of the Hi 
(21 -cm line) interstellar gas distribution on a 40^ x 40^ field around 
Geminga based on VLA and Effelsberg radio telescope data. Based 
on the line-free channels of the 1 .4 GHz cube, they produced a ra- 
dio continuum image with 24^^ angular resolution, confirming the 



^ See http ://www. aips .nrao . edu/ ' 

^ For more details, see the VLA observation status summary at 
|http ://w w w. via, nrao . edu/ astro/ guides/vlas/ current/ 1 



^ We checked our capability to distinguish point-like from slightly elon- 
gated sources, such as the detected feature, by performing detailed simula- 
tions in which we artificially injected point sources in the raw VLA data. 



Continuum radio emission associated with Geminga L3 



nondetection of any source at this frequency, either point- hke or 
extended, down to a noise level of 0.14 mJy beam~^ (la). 

Hence, the presently available data do not allow one to con- 
strain the thermal or non-thermal nature of the emission, but, if the 
radio source is steady, we can set an upper limit on the spectral in- 
de^lflof af l §Hz < 0.1 (this value has been calculated by consider- 
ing that at 1.4 GHz the total flux density is below 3a i a ghz = 0.42 
mJy). We note however that if the radio source is variable, this limit 
is unreliable, since the reported observations are not simultaneous. 

We computed the coordinates of the Geminga pulsar at the 
epoch of the reported VLA observation by projecting the ac- 
curate optical absolute position by Caraveo et al. (1998) using 
the proper motion measurement by Ipahertv et al.l (I2007D : RA = 
06^33"'54!24, Dec. = +17°46'13'.'9 with errors <0'.'l. The pul- 
sar is within the reported diffuse feature, marginally displaced by 
2.7 db 1.8 arcsec from the radio peak emission and in association 
with a radio brightness of 0.26 mJy beam~^ (^-^80% of the peak 
flux; Fig.lJJ.The probability of random association between the ra- 
dio feature and the Geminga PWN region (approximately FWHM 
beam size divided by the field of view) is <10~^. We examined 
radio, optical and near-infrared catalogs looking for possible back- 
ground objects of different nature (e.g. an AGN jet), but no coun- 
terpart compatible with the observed structure was found. 

In order to compare the radio and X-ray morphologies of 
the structures associated to Geminga, we retrieved a public Chan- 
dra ACIS-S imaging exposure o f the field p erformed clos e to the 
VLA observation (Obs. ID 4674: iDe Luca et a l. 2006: Pa vlov et al.l 
l2006h . The ^20-ks-long observation was carried out on 2004 
February 07 (168 days before the VLA observation) in 'Timed Ex- 
posure' mode using a 1/8 subarray. Geminga was positioned in the 
back-illuminated ACIS-S3 chip, sensitive to photons between 0.2 
and 10 keV. The data were processed using the Chandra Interac- 
tive Analysis of Observation software (CIAO, version 4.2) and we 
employed the CALDB 4.3 calibration files. Standard screening cri- 
were applied to generate a 0.5-8 keV image of the field of 
Geminga. Using the Chandra aspect tooljf] we found a small offset 
in the astrometry (ARA = O'.'Ol, ADec. = O'.'l). After fixing the 
offset, in order to assess the absolute astrometry of the Chandra 
data set, we computed the position of Geminga with the WAVEDE- 
TECT task. We localized it at '-^O^.^S from the position expected on 
the basis of the source absolute optical position and proper motion, 
well within the typical ACIS-S localization accuracyO The X-ray 
contours obtained from the 2004 Chandra observati on and related 
to t he pulsar and to the axial PWN tail (see also iDe Luca et al.l 
I2OO6 : Pavlov et al. 2006^) are shown in Fig. |2] The deconvolved 
length of the radio feature (^10^0 is smaller than the whole time- 
integrated (2004 and 2007 Chandra combined data) X-ray tail, but 
larger than the substructures (typically ' ^5^^-long 'blobs') detected 
by Chandra in individual observations jPavlov et al.ll2010l) . 



3 DISCUSSION 

A VLA observation at 4.8 GHz unveiled a ^10^^ -long radio feature 
with flux of ^-^0.4 mJy associated with the Geminga pulsar and its 
PWN. Significant contamination from Geminga persistent pulsed 
radio emission is very unlikely: given the tight upper limits on 
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Figure 2. Detail of the radio feature shown in Fig.l. Black cross marks 
Geminga pulsar position at the epoch of VLA observation (July 2004) and 
the dashed arrow indicates its proper-motion direction. Green contour levels 
represent the X-ray emission (pulsar and PWN) detected by Chandra in 
Febru ary 2004 (see Section |2). Most of X-ray axial features jPavlov et al.l 
l20ld) are detected within ~10^^ ahead and behind the radio peak emission. 

pulsed flux density at low frequency (e.g. < 0.05 mJy at 0.4 GHz; 
iBurderi et al.lll999h . the contribution of pulsed emission to the to- 
tal observed flux at 4.8 GHz would be <0.1% for the mean spec- 
tral index of radio pulsars «a>= 1.8 =b 0.2: lMaron et al.ll2000h 
and could at most reach ^15% even in the extreme case of a flat 
spectrum. The occurrence of transient pulsed (or continuum) radio 
emission in isolated NSs in not unheard-of, but it involves at our 
present knowledge pulsar classes showing a totally different phe- 
nomenology from that o f the Geminga pulsar, l ike magnetars and 
rotating radio transients (iRea & EspositolEoill : iMcLaughlin et al.l 
l2006h . Furthermore, the emission models mentioned above, as well 
as any other kinds of magnetospheric emission originating very 
close to the NS surface (and hence appearing as point-like at the 
spatial resolution of the available observations), are disfavored by 
the evidence for a slightly elongated feature, together with the oc- 
currence of a marginal offset between the NS position and the radio 
emission centroid (~ 2^.^ 7). 

The latter two features and the tendency of the radio emis- 
sion to distribute along the direction of the axial X-r ay tail, natu- 
rally l ead to an association with the Geminga PWN. IPavlov et ahl 
(I2OIOI and references therein) observed that the axial tail may re- 
sult from the shocked pulsar wind collimated by the ram pres- 
sure (i.e. a bow-shock) or from a jet emanating along the pulsar's 
spin axis and perhaps aligned with the direction of motion. In the 
frame of the standard bow-shock models, radio PWNe are gener- 
ally interpreted as arising from cooling of X-ray-emitting leptons 
jKennel & Coronitill 19841) . The smaller size of the radio nebula of 
Geminga (with respect to the whole, time-integrated, ^50^^ long X- 
ray tail) apparently contradicts this hypothesis, unless we associate 
the observed radio feature to the variable streaming 'blobs' of ^^5'^ 
diame ter, possibly consti tuting the X-ray axial tail as observed in 
2004 jPavlov et al.ll2oTQh R Furthermore, synchrotron models ex- 
plaining the X-ray tails as pulsar wind shocked by the ram pressure 



^ S(v) oc where a is the spectral index. 

^ See http://asc.harvard.edu/ciao/threads/index.html. 

^ See http://cxc.harvard.edu/cal/ASPECT/fix_offset/fix_offset.cgi. 

9 See |ht^://cxc.harvard.edu/cal/ASPECT/celmon| 



Note that since presently available Chandra and VLA observations are 
not simultaneous and the X-ray axial tail is partially hindered by the bright 
NS, we cannot in any case unequivocally link the observed radio feature to 
a specific substructure of the tail. 
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due to the NS's supersonic moti on typically require r elatively low 
magn etic fields (-10-100 uG: ICaraveo et all l2003l : IPavlov et all 
l2006h . With such low magnetic fields, the radio emission arising 
from the cooling of the X-ray emitting electrons, would be placed 
well outside the VLA field of view due to the large proper mo- 
tion of Geminga: e.g. more than 0.5 deg off the position of the 
X-ray emission for a magnetic field of 100 /iG|^ In summary, 
ram-pressure collimated PWN models can be applied to our case 
only if additional processes, such as adiabatic expansion from hy- 
pothetical over-pressurized substructures of the tail (not observable 
at the present stage) provide a cooling power much higher than the 
synchrotron one and/or if radio emitting electrons are directly ac- 
celerated at the PWN shock, i.e. they are not originated from the 
cooling of high-energy electrons as in standard PWN models. In 
this respect, we note that the radio "wisps" observed in the Crab 
PWN are produced by the same process as the higher-energy ones, 
although the mechanism by w hich this particle spectrum is gener- 
ated is still not known (see e.g. Ib ietenholzll 200^ . 

Alternatively, the radio tail of Geminga could in principle 
be ascribed to a (magnetically-collimated) jet, in which the ra- 
dio emission could directly result from radio-emitting electrons 
freshly-injected along the spin-axis at the jet formation site, or 
from internal shocks in which the accelerated X-ray-emitting lep- 
tons cool down to radio energies on the observed short space and 
time scales ( see e.g. lBenfordll984l : rKomissarov & Lvubarskvl2004l : 
lMeieill2003l for possible models of magneto-hydrodynamic jet pro- 
duced by isolated pulsars). 

In the latter case, magnetic fields >0.05 G are required in or- 
der to avoid the occurrence of an additional (i.e other than syn- 
chrotron) cooling process. As already noticed for the case of the 
ram-pressure collimated PWN models, there is no evidence of such 
high fields in other PWNe. However, in the first place, large mag- 
netic fields due to the NS dipole are in principle available at the 
pulsar light cylinder {B\c ^ 10^ G for Geminga) and in its sur- 
roundings. In the second place, an observational bias may also be 
at work here. In fact, despite the grow ing number of jet-like struc - 
tures around NSs resolved in X-rays (iKargaltsev & Pavlovll2008h . 
no 'genuine' magneto-h ydrodynamic radio jet around pulsar was 
detected so far (see e.g. ICohen et aDll983l for VLA mapping of 
numerous pulsar fields), possibly excepting PSRJ202 1-1-4026, an- 
other radio-quiet gamma-ray pulsar, located at the edge of the su- 
pernova remnant G 78.2-1-2.1 (iXrepl et al.ll2010l) . Thus, it cannot be 
excluded that pulsar radio jets might be observable only in nearby 
and radio-quiet pulsars, not hindering their close environment by 
bright magnetospheric radio emission and having magnetic fields 
much stronger than those related to other PWN structures placed 
at larger distances from the NS surface. Interestingly enough, this 
hypothesis may be checked by a series of deep radio interferomet- 
ric exposures. Scaling our result for Geminga, this survey could 
unveil jet-like structures in few other near (<1 kpc) and energetic 
(^rot ^ 10^^ erg s~^) radio-qui et pulsars, some o f them recently 
discovered by the Fermi satellite (lAbdo et al.ll2009h . 

Finally, we note that a further possible explanation of the 
Geminga radio tail is the leakage of radio-emitting electrons (with 
Lorentz factor 7 ^ 10-100) directly from the open field lines of the 
magnetosphere. The maximum current through the mag netospheric 
accelerator (i.e., the Goldreich- Julian current; see e.g. IWang et aP 



See also e.s. I Wane et al.ll200ll and lOaensler et al.ll2004l for examples 
of parameters of typical bow- shock PWNe simultaneously detected both in 
X-rays and radio. 



1 19981 and references therein) and then the maximum particle flux 
escaping from Geminga open field lines can be estimated from 
N ^ Q'^BsurfR^/2ec ^ 5 x 10^^ s"^ (where Q is the pulsar 
angular velocity, and the NS's radius R is assumed to be 10 km). 
The expected synchrotron power at 4.8 GHz, for a random distri- 
bution of pitch angles, can be roughly obtained from: 

^rTdio^ — ^<^T^^UBcNeTsynch ^ 10^^76 erg 

where Tsya=mecl {4:(Tt^U b) is the synchrotron cooling time of the 
electrons. Us — j (Svr) is the magnetic field energy density and 
e is the fraction (w.r.t. iV) of radio emitting particles. The observed 
luminosity of 10^^ erg s~^ at 4.8 GHz would then be matched 
for 76 10, implying a magnetic field oi B ^ 20e^ G in the 
area of the emission. If this magnetic field is associated to the NS 
dipole, the emitting region would then be placed at a distance of 
^em ^ 5 X 10^e~^/^ cm from the NS surface. In summary, for 
e = 0.5 (in turn implying Dem of the order of ten light cylin- 
der radii, 7 20, 5 ^ 5 G, and Tsyn ^ 0.05 yrs), this model 
would provide a cooling and possibly expanding diffuse radio fea- 
ture protruding from the Geminga pulsar. If this model is correct, 
two geometrical constraints should be simultaneously satisfied by 
the Geminga pulsar: the radio pulsar beam not intersecting the line 
of sight and the leaking electrons producing a relatively narrow 
cone of radio emission roughly aligned with the projection of the 
pulsar proper motion. This would suggest that Geminga is a nearly 
aligned rotator, with the spin axis almost aligned with the proper 
motion and with the velocity vector close to the plane of the sky. 



4 CONCLUSIONS 

We have reported on the discovery of a radio emitting structure 
with a flux of ^-^0.4 mJy at 4.8 GHz, associated with the Geminga 
pulsar and roughly aligned with the direction of the neutron star 
proper motion. The source appears slightly extended (^^ 10^^ -long) 
but further deeper observations are needed to confirm the elonga- 
tion unambiguously. The nature of this emission is unconstrained 
by present data, but a combination of dedicated radio interferomet- 
ric and X-ray observation£3 is expected to shed light on the under- 
lying physical process and likely unveil similar radio tails in other 
radio-quiet pulsars. 
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In particular, bulk flow velocity measurements in the radio band could 
assess the indication of motion/variability of X-ray blobs along the axial 
tail, suggestive of a jet-like inter pretation, obtained by comparing Chandra 
observations at different epochs dPavlov et al ]|201Q|) . 
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